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Abstract 

It is explained how to provide self-adjoint operators having scattering states forming a mul- 
tiplicity one continuum and bound states whose corresponding eigenvalues have an asymptotic 
density equivalent to the one of the zeros of the Riemann zeta function. It is shown how this 
can be put into an integro-differential form of a type recently considered by Sierra. 



1 Introduction 

The present publication was motivated by a recent investigation of Sierra [23] (also [2H [21]). It 
will be explained how the differential system and isometric expansion we obtained in [7J, and the 
associated self-adjoint operators, provide, when moved back to the "rr-picture" a particular example 
of a type of differential-integral equation considered in [531 section IV] . 

In [5] we determined "explicitely" certain objects depending on a parameter a > and associated 
to the cosine (or sine) kernel. In [7] we obtained the differential system (which turned out to be 
of the Dirac type with coefficients involving some Fredholm determinants of the cosine (or sine) 
kernel) obeyed by the objects from [6] as functions of log(a). We mentioned the associated isometric 
spectral expansion and also explained how this investigation had led to the realization of the Fourier 
transform as a scattering, an objective which had resulted from our definition and study of the 
"log |x| +log|p|" operator ([!]). 

In fTU] we gave a detailed exposition of these results. Some general aspects could be predicted in 
advance as they involve some Hilbert spaces of entire functions in the sense of [3J . The specific spaces 
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had actually already be defined in [2j [20], among spaces associated to general Hankel transforms. 
The results from [5J [5D] regarded mainly the spaces associated to the Bessel function Jo. In [H] 
we obtained for the cosine and sine kernels the "explicit" form taken by some functions whose 
existence was a consequence of the general theory of [3]: in particular certain entire functions A a (s) 
and B a (s) having all their zeros p on a line, and such that the quotients ^ai£j (resp. ) give 

orthogonal bases of some Hilbert space K a . It was then shown in [7] that these objects obey 
differential equations whose coefficients arc expressed in terms of certain Fredholm determinants. 
This observation also applied to the Jo kernel and we gave the detailed exposition of these results in 
[TO] , where the chapters V and VI with minor adaptations apply to all Hankel kernels, in particular 
to the cosine and sine kernels. A further advance had been realized in the theory of the Jo kernel 
and associated spaces whose nature was revealed to be special relativistic [H] . Other Hilbert spaces 
of entire functions, directly associated with the Ricmann zeta function and which had not been 
considered before, were defined by the author in [5] and shortly thereafter understood to be in 
connection with the "co-Poisson" formula. 

In [23j , Sierra considers integro-differential operators which are singular rank two perturbations 
of + |) on L 2 (0, +oo; dx). We will explain how the objects, differential equations, and 

spectral expansion from [7] [TU], when seen in the u x picture" lead to the kind of equations 
considered by Sierra. As was mentioned in a number of our other papers (such as [8]) the theory 
leads to "bound states" associated to zeros on the critical line which have an asymptotic density 
equivalent to the one applying to the zeros of the Riemann zeta function. These bound states 
are created on an interval (log(a), +oo) by a Dirichlet condition at log(a) and a potential function 
which has an exponential increase at +oo. The Dirac system seen on (— oo, log(a)) leads to a 
purely continuous spectrum, because the potential vanishes exponentially quickly at — oo. Hence 
we do have here the combination of scattering and bound states mentioned in section IV of [23] 
as a possible property of the integro-differential equations considered there. It is not clear if the 
particular example which is considered in section VI of [23] truly leads to bound states, or rather 
only to resonances. On the other hand, the Dirac and Schrodinger equations from [7] do have, as 
we mentioned in our publications, this property of leading to "quantum zeros" lying on the critical 
line with the expected density (a general density result has been given in [S], see Theorem 7.7 and 
Remark 19). 

We shall not enter here into other topics considered by Sierra, such as obtaining operators leading 
exactly to the zeros of Riemann. 
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2 Some notations 



For the sake of facilitating comparisons with [23] we shall often write s = \ + iE when using complex 
numbers although the letter E is also used with a completely different meaning. Also, our scalar 
products (■ | ■) will be conjugate linear in the left entry (bra) and complex linear in the right entry 
(ket). As we refer to somewhat lengthy developments from our previous papers we did not modify 
our notations too much, in particular with respect to the use of the letters A and B. We will not 
here attempt to analyse our final equations in the spirit of [231 IV] but only want to show how they 
emerge from our framework, so the incompatibility of notations will remain virtual. 

The Mellin transform is given by the formula /(| + iE) = / °° f(x) x ^ dx — / °° f(x)x~ s dx 
(and not with +iE) because this represents the scalar product (ip%\f), of / with the generalized 
eigenvector if>%(x) = ^= of the "free" operator H = + i), Ha(ip%) = Etjjg. The inverse 

Mellin transform is given by the formula f(x) = ^ f(\+iE)^ dE = ± Jsr( s )=i /(«)^ s_1 \ds\. 

Regarding the Fourier transform, it will be used in the form of the cosine transform C(f)(y) = 
/ +O ° 2 cos(2TTxy)f(x)dx, for functions or distributions supported in [0, +oo). The cosine kernel 
2cos(27ray) has only the two eigenvalues +1 and —1 on (0, oo). When restricted to a finite interval 
(0, a) it has a discrete spectrum 1 > |Ai(a)| > | A2 (a) | > ••■ which accumulates at zero. This 
restriction is conjugated to, hence has the same spectrum as, the kernel C a (x,y) — 2acos(2ira 2 xy) 
acting on the interval (0, 1). We will make use of the Fredholm determinants det(l ± C a ). Some of 
their properties (in particular the asymptotic as a — > +00) can be derived from knowledge of the 
corresponding facts for det(l — C%), As the square C\ is the Dirichlet kernel 8m ^°J^~ 1/ ^ acting 
on the even subspace of L 2 (— 1, +1; dx), the question is reduced to the finite Dirichlet kernel, whose 
properties have been extensively studied in the literature in connection with random matrices: see 
[P7] and the further references included in our papers [jj [10] . 



3 A support condition and associated Hilbert spaces 

To each square integrable function / on (0, +00), one can associate its Mellin transform f(s) = 
Jo°° f{x)x~ s dx on the critical line 3?(s) = \. This Mellin-Plancherel transform is a unitary iden- 
tification of L 2 (0, +00; dx) with i 2 (5R(s) = i; ^-). The closed subspace L 2 (l, +00; dx) (resp. 
L 2 (0, 1; dx)) becomes identified with (boundary values of the functions of) the Hardy space of 
the half-plane 5R(s) > \ (resp. 5i(s) < |). Let us recall how the cosine transform g(y) = C(f)(y) = 
J +o ° 2 cos(2w xy)f(x)dx is represented on the Mellin side. The formula is: 

C{f)(s)=x(s)f(l-s) where X ( S ) = ^~*^!^ W 
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This is the same x(s)p of modulus one on the critical line, which appears in the functional equation 
C( s ) — x( s )C(l — s ) °f the Riemann zeta function. Indeed this functional equation can be seen to 
be a manifestation of the self-invariance under the cosine transform of the distribution ^2 n>1 ~ 
n) — l x> o(x), whose Mellin transform is nothing else but £(s) itself. 

Of course £(s) is not square integrable on the critical line but the functions s i— > are, with p 
a zero of the zeta function. Let us consider the vectors Cp,fc( s ) — (s- S p) k wnere P is a non-trivial zero 
of zeta, and k is an integer between 1 and the multiplicity of p. Taking inverse Mellin transforms 
we also tacitly consider the £ p ^ to be vectors in L 2 (0, +00; dx). The following holds ([SJ): 

Theorem 1. The vectors are linearly independent in the sense that none is in the closure of 
the linear combinations of the others. Furthermore the smallest closed subspace of L 2 (0, +00; dx) 
containing these vectors consists exactly of those functions f in L? which are constant on (0, 1) and 
whose cosine transforms are also constant on (0, 1). 

We mentioned this theorem as motivation to study the space K% of functions / in L 2 which 
vanish identically on (0, 1) and whose cosine transform also vanish identically on that same interval 
(the condition of vanishing is slightly easier than the one of being constant). More generally we 
can consider the space K a where the interval (0, 1) is replaced with the interval (0,a) (we would 
have been led directly to these spaces if we had considered the functions Ll fj^ > where x is an even 
Dirichlet character): 

Definition 2. We let, for each a > 0, K a be the subspace of L 2 (0, +00; dx) consisting of functions 
which vanish identically on (0, a) and whose cosine transforms also vanish identically on that same 
interval. 

The subspace K a is closed and its orthogonal complement is the (non-orthogonal, but direct and 
closed) sum L 2 (0, a; dx) + C(L 2 (0, a; dx)). 

Proposition 3. If f € K a then f(s) is an entire function which has trivial zeros at s = 0. —2. —4, 
-6, . . . 

Indeed f(s) as Mellin transform of / which is supported in [a, +00) is analytic in the half-plane 
3t(s) > h, and its boundary values on the line are the Mellin-Plancherel transform of / (taken a 
priori in the L 2 -sense). The same applies to g = C(f) and on the critical line one has the almost 
everywhere identity g(s) — x( s )f0- ~ s )- So the function k(s) — x(s)g(l — s) is analytic in the 
half-plane 5R(s) < |, and its boundary values from the left coincides with the boundary values of f 
from the right. From general principles of complex analysis ( "edge-of-the- wedge" theorems) the two 

Met us observe in passing x( s )x(l — s) = I and using formally f(x) = 8(x — 1) that x{ s ) = /o°° ^ cos(2nx)x~ s dx. 
This is valid as a semi-convergent integral, and also as an identity of tempered distributions in log(x) and ^s(s) 
(«(») = i). 
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analytic functions "glue" to give an entire function. It would be easy to give in this one-dimensional 
case more details for this glueing, but anyhow various simple elementary constructive proofs of the 
Proposition are available, such as the one given in [5 The trivial zeros are a corollary of the 
identity C(f)(s) = x( s )/(l — s )j because /(l — s) must have the zeros to compensate the poles of 
x(s) at 1, 3, 5, . . . 

The spaces K a , mainly as spaces of the entire functions J-(s) — 7r _3 r(|)/(s) (but seen as 
functions of the variable z such that s = i — iz), had appeared earlier in the literature: in de Branges 
investigation [2] where they are studied as a special instance of spaces verifying the three axiomatic 
propertied whose general consequences are the subject of the book [3]. The well-known Paley- 
Wiener spaces (of Fourier transforms of L 2 functions of compact support) are another, simpler, 
instance of the general theory of [3]. A particularly interesting (and challenging) aspect of the 
spaces K a is that the associated entire functions are not of finite exponential type. The reader 
interested generally speaking in [3] is also referred to [H] and to [13] and [19], especially for the 
rather intimate connection with standard aspects of spectral expansions associated to differential 
equations of the Sturm-Liouville or Schrodinger types. 

In the process in [6] and [7] of studying the spaces K a directly we have in particular recovered 
by other means the consequences predicted by the general theory of [3], in an "explicit" form. In 
this section and the next we briefly review some of these results. 

The evaluation maps / i— ► /(s) are continuous linear forms on K a , they thus correspond to 
certain vectors in K a indexed by the complex numbers. The space K a is "elucidated" once the 
mutual scalar products between these evaluators are computed. Let Z z be the entire function such 
that T(z) = (Z z \3-) for all completed Mellin transforms J- of elements of K a (we may consider 
Z z either as an element of K a C L 2 (0, oo; dx) or in its incarnation as an entire function Z z (s); 
the dependency on z is anti-analytic). Then, from J3] it is known that there must exist an entire 
function £ a (not quite but almost unique) such thato 

z z (w) = (z w \z z ) = W£ a W-gg^)g a (i-^ (2) 

Such an £ a once found also characterizes the space of the entire functions !F(s) as follows: the 
ratio y~ must belong to the Hardy space of the half-plane 5ft(s) > ^ and also the ratio f- where 
Q(s) — JF(1 — s). The formula ||/||2 = II^IIh 2 holds, that is, somewhat surprisingly: 

-1 / |/(i + iE)f dE= — [ \H- + ^)| 2 - - I ffl-^' 2 dE (3) 

2Wk 1 2 U 2tt 7 M 1 4 n \£ a {\+iE)\* 2nJ R \ Ea (l +i £)| 2 

2 in [5] we use f(x)x s ~ 1 dx whereas systematically in later references we always use Jq°° f(x)x~ a dx. Hence the 
s from [5] is the 1 — s from here. 

3 not needed here, see [3], and for the K a spaces [2l I20| and [6l l8l [To] ■ 

4 we replace here the horizontal line always used in [3] by the critical line = ^. Hence formula Q does not 

look as in [3]. 
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We have introduced here the notation E a (s) = £ a (s)/7r ^r(|) and indeed this proves useful when 
studying the spaces K a . 

The reader is referred to [3] and P3] for the general proof (which is short) of existence of such 
"^-functions" when the axioms of [2] are verified. From the evaluator identity one deduces that 
|£ a (s)| 2 > |£ a (l — s)| 2 when 5R(s) > |, so all zeros of £ a {s) verify 5R(s) < h. A zero on the 
symmetry line would mean that all functions .F(s) vanish there. It can be shown that this is not 
the case, so the zeros of £ a (s) in fact verify Sft(s) < \. 

The author has contributed the "explicit" determination of a suitable function £ a |6J. It is 
realized as n~^T(^)E a (s) where E a {x) is in fact a distribution supported on [a, +oo) and whose 
cosine transform is also supported on [a, +oo). To obtain E a , we studied the vectors X s (x) which are 
such that f(x)X s (x) dx = f(s) for / <G K a . For 5ft(s) > \ the vector X s is simply the orthogonal 
projection to K a of x~ s l x>a (x). We determined directly the integrals f°° X z (x)X w (x) dx. A 
completely detailed exposition with all proofs has been given in our manuscript [10 . As explained 
there, the reader will only need to replace all occurences of the kernel Jo(2y/xy) by the kernel 
2cos(27rxy), and also to replace all occurences of T(s) by tt~^T(^). The chapters of [TD] which are 
related to the present publication are the fifth and the sixth. The simplifications allowed by the 
kernel Jo(2^/xy) are made use of only in other chapters of [TP] . 

The formula obtained in [6] is: 
Theorem 4. Let <fr^ be the even entire functions solving: 

( t ) i( x )+ / 2 cos(27rxy )</>+(?/) dy = 2cos(27raa;) (4a) 
Jo 

r 

4> a ( x ) ~ f 2cos(27rxy)0 Q (y) dy — 2cos(2nax) (4b) 



o 



The function E a (s) given as 

E a {s) = ^ [a- s + ~ J™(4>i(x) - r a {x))x- s dx^j (5) 



is an entire function and is such that for all z,u)GC 



oo 



= A;;, = / XWX^dx = ^^^^^1 (6) 



a 



Equation © looks different from (0) because we dropped Gamma factors and also because we 
are not using the Hilbert scalar product to avoid anti-analyticity. Equations (|4"aj) . (|4b[) are first 
considered in L 2 (0,a;dx). They have solutions because ±1 are not in the spectrum of the finite 



5 Lagarias has considered in 1151 an entire function which has the properties of an ^-function a la de Branges if 
and only if the Riemann hypothesis holds. Lagarias E function would vanish at the locations on the critical line of 
the non-simple zeros of zeta, if some exist. 
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cosine kernels. They imply that 4>^{x) and <p a (x) are even entire functions, in particular they have 
meaning for x > a. The integral for E a (s) is absolutely convergent for 3?(s) > 0. 

The function £ a (s) is real for s real. Writing £ a (s) = A a {s)—iB a (s) on the critical line, with A a (s) 



and B a (s) real there, we thus have for s on the critical line: 2A a (s) — £ a (s)+£ a (s) — £ a (s) + £ a (s) = 
£a{s) + £ a iX — s )- The last equation, for general s, defines A a as an entire function, which is even 
under s i— > 1 — s, real for s real, and real for s on the critical line. Also 2B a (s) — i(£ a (s) — £ a (l — s)) 
is an entire function, which is odd under s i— > 1 — s, imaginary for s real, and real for s on the critical 
line. We shall write A a {s) = 7T~ir(|)A a (s), B a (s) — 7C~%T(^)B a (s), so that A a (x) (resp. B a (x)) is 
a cosine invariant (resp. skew-invariant) distribution supported in [a, +oo). With 5 a (x) = 6(x — a) 
one has indeed: 



Aa{x) = ^(8 a (x 
-iB a {x) = (S a (x 
The evaluator identity $2$ can also be written as 

\2 u: j Z 

In particular, for z — w one obtains 

IIZJI 2 



^(2 ; ) 1 a;>a(a;)) 
C( a; ) 1 x>a(a;)) 



A a {z)B a {w) - B a {z)A a {w) 



i(z + u> — l) 



MA a (z)B a (z)) 



\Zi +iE \\ 2 



.1 



1 



-2A(^ + iE) iB' a {- + iE) + 2B a (- + iE) iA' a {- + iE) 



(7a) 
(7b) 

(8) 

(9) 
(10) 



Derivatives in (jTUJ) are taken with respect to s = i +iE (not E). Equation ^ implies that the entire 
functions A a and B a have all their zeros on the critical line. It is indeed a fact that all vectors Z z 
are non-zero (it is easy to see that "T(z) = for all JF" is impossible). From we deduce that the 



ratio 



A a {s) 



is of positive imaginary part when 9?(s) > |. It is elementary, as 



is meromorphic, 



that this implies that its poles and zeros (all on the critical line) are simple. Next, if A a had a 
multiple zero p, then B a would have to vanish there, but then taking a limit in ^ we would obtain 
\\Z p \\ = 0, which we have already mentioned is false. So the zeros of A a and B a are simple; using 
()10p one argues that the respective zeros interlace. These arguments apply in the general situation 
of 0. A simplification which is valid here and that we have made use of is that the evaluators Z z 
are non-zero also when z is on the symmetry line. 

From ([5]) two vectors Z p and Z p * associated to distinct zeros of A a are always orthogonal. 



Furthermore from 
Pn = 1 - Pn 



if p n = 2 + iE n is a zero of A a then Z Pn (s) — 2iB a (p n 



So the vectors are mutually orthogonal. From (flU|) : 



S-Pn 



A a (s) 


' "I 


A a (s) 


2 \ds\ 


1 d 


Aa ( 2 


V iE) 


S - Pn 


K a «(s)=i 


S - Pn 


2vr 


~ 2 dE 




V iE) 



(we used 



(11) 



7 



We can write this also as 



A a (s) 


1 


d 




S - Pn 




dE 


E=E 



A a a + iE) 



E =E n B a a+iE) 2 idE 



Ml+iE) 

E=E n Ea{\+iE) 



(12) 



The vectors corresponding to the zeros of B a provide another orthogonal system. These vectors 
corresponding to the zeros of A a (or B a ) span all of K a . This orthogonal basis theorem is a general 
fact from [2j (under circumstances which are verified here). 

In the next sections we prove this orthogonal basis theorem from the standard material provided 
by the Weyl-Stone-Titchmarsh-Kodaira theory of second order differential equations ([TU HB]), as 
an outcome of studying the deformation with respect to a > of the spaces K a ([7]). Anyhow, 
obviously one of the main origins of the axioms of [3] is precisely the behavior of the objects from 
the Weyl-Stieltjes spectral theory of differential or difference equations (see Remling [T!5] for the 
spectral problems of Schrodinger equations from the point of view of these axioms). 



4 Dirac system and associated isometric expansion 

Definition 5. We let C a be the kernel 2 cos(2nxy) acting on L 2 (0, a; dx). And we define 

with u = log(a) 



, » d . det(l + C a ) 
H(u) = a— log ■ 



da ° det(l - C a ) 

Definition 6. We let A(u, s) = A a (s) and B(u, s) = B a (s) for a = exp(u), — oo < u < oo. 
Theorem 7 ( 7 ). The following differential system holds, where s = ^ + iE, E G C: 

(— + fi(u))A(u, s) = -EB{u, s) 
au 



{— - n(u))B(u, s) = +EA(u, a) 
du 



Equivalently: 





' 1" 


d 





u(u) 


) 


A(u, s) 


= E 


A(u, s) 


( 


-1 


du 


ji{u) 







B(u,s) 


B{u,s) 



(13) 



(14a) 
(14b) 

(15) 



Combining the differential system (flS]) and the evaluator formula (|8]) one obtains the following 
identity 



POO 

{Zl\Zl) = / 2A{u, z) 2A{u, w) + 2B(u, z) 2B(u, 

Jloe-(a) 



w) 



du 



(16) 



og(o) 

and in particular we sec that for each z G C the two-component vector 
at +oo. 



A{u, z) 
B(u, z) 



is square- integrable 



Let us recall that according to a general theorem [T^l §13, Thm 7.1] for Dirac systems with 
continuous coefficients we are automatically in the limit-point case at infinity. Being in the limit 



8 



point case at +00 means that for each E £ K there is exactly one (up to a scalar multiple) non- 
vanishing solution which is square-integrable at +00. Thus we have identified this solution, and 
we have seen that it also exists for E real (the solutions which are square-integrable at — 00 are 
also known [7]; we will mention them in a later section). There can't be other linearly independent 
square-integrable solutions as this would characterize the limit-circle case, which never happens for 
Dirac type systems with continuous coefficients!^ 



Let us consider the Dirac-type differential operator 



" 1" 


d 





fi(u) 


-1 


du 


ji{u) 






acting on pairs 



01 { u ) 

; ~°° < u < +00. As we are in the limit-point case at infinity (plus or minus), from Hermann 
Weyl's theory (adapted to first-order two-component systems) we know that the symmetric operator 
defined by this differential operator with domain given by the functions of class C 1 and with compact 
support is essentially self-adjoint ( [161 [TT1 [T5] ) . In Theorem 15 of [TO] we determined the associated 
isometric spectral expansion and showed that this self-adjoint operator Hq is — i{x-4- + \) acting on 
L 2 (0, 00; dx) or equivalently the operator of multiplication by E on the critical line: 



Theorem 8 (10 , Theorem 15). There are unitary identifications between L 2 (0,oo;dx), L (5i(s) = 
1. kM \ ana i £2qjj £2. du-^ given by the following formulas (to be understood as suitable 



2 ' 2;r|7r-7r(§ )| 2 

L 2 -limits, see IWj), where always s = \ + iE, E 6 



a(u) 



\f{x)\ 2 dx 
Under these identifications: 



f(x)x 

<X<OQ 

2A(u, s) 
2B{u,s) 



' dx 
dE 



a{u) 2_A(u, s) 



27T |7T-2r(§)| 2 
nil 

(3(u) 2B(u,s) Y 



rap 



dE 



27r|7r-4r(i 



(\a(u)\ 2 + \P(u)\ 2 )^ 



(17) 
(18) 
(19) 
(20) 



• The orthogonal projection of f(x) to K ao corresponds to 



a(u) 




a{u) 


J3(u)_ 


1 — > 


JS(u) 



l«>log(o )W> 



the Cosine Transform f 1— > C{f) corresponds to [p] 1— > [—p] (and F(s) 1— > ^"(1 — s)), 

the self-adjoint operator —i(x-^ + ^) on L 2 (0,+oo;dx), which multiplies by E with maximal 
domain on the critical line, corresponds to the unique self-adjoint completion of the differential 
operator on the u-line: 



Hn = 



4." 

du 









—f 

du 










(21) 



with initial domain the functions of class C 1 ( or even C°° ) with compact support. 



6 in |10| we give an eigensolution for E = which we prove to be not square-integrable at infinity. This explicit 
argument confirms that we are not in the limit-circle case. 
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This Theorem is an elaboration deduced mainly from the evaluator identity (I16|) . Let us also 
mention that the integrals in (|18[) are absolutely convergent when / is in the domain of Hq, that 
is, when s ■ f(s) is square integrable on the critical line. The functions a(u) and (3(u) are then 
continuous; actually we know that they are absolutely continuous as they must be in the domain of 
Ho in the u-picture, so their derivatives must be locally square integrable. 

The operator xp = —ix-£^ : or \{xp + px) — —i{x4z + \) was considered "semi-classically" 
by Berry and Keating :T in a discussion where it was accompanied with the conditions \x\ > l x , 
\p\ > l p (the operator —i(x-4- + |) with conditions \x\, \p\ < A arose in the work of Connes [H]). 
Certainly it is natural to consider that the conditions |x| > l x , \p\ > lp, treated "semi-classically" 
in PQ, could at the "quantum" (that is, operator-theoretical level) be related with the condition 
of support which defines the spaces K a (and if we were considering functions odd in x, we would 
have the analogues of the spaces K a associated with the sine kernel). Recently, German Sierra in 
[23] (see also [2UE2]) has generalized the conditions |x| > l x , \p\ > l p through the consideration of 
differential-integral equations obtained from the differential operator xp by the addition of (singular) 
finite rank perturbations. Sierra presents general aspects of the corresponding eigenfunctions. This 
leads to scattering states (deformations of the "free" states x lE x~^), and possibly also to bound 
states in some fine-tuned situations. In the next section we see how we have natural bound states in 
K a and then we will show that in the rr-picture these states are eigensolutions in the form considered 
in [33] Section IV]. The scattering states will span the perpendicular complement of K a in L 2 . 



5 Bound states 



We consider the differential system ([15]) on the half-interval (uq, +oo), uo = log(do) with boundary 
condition a(u ) — 0. From general theorems [TSJ [TBI E] we know that this leads to a self-adjoint 
operator and associated isometric expansion (with discrete spectrum from the behavior of fi(u) at 
+oo, but we will confirm this directly). To obtain the spectrum we apply Hermann Weyl's method 
( [111 §9], [IS1 §3]). For each s E C we let ip(u,s) be the unique solution of the system (fT5)) for 
the eigenvalue E, s = \ + iE, and with the initial condition i/j( u o, s ) — [J] an d ^ 4*{ u i s ) be the 
unique solution for the eigenvalue E with the initial condition </>(ito, s) = [\]. Let, for ^s(E) > 0, the 
quantity m(E) be chosen such that ifi(u, s) — m(E)<j>{u, s) is the unique linear combination which 
has the property of being square- integrable at +oo. 

To check up on signs, let us see what would happen if the potential /Lt(u) were identically vanish- 



ing. Then we would have ^(u, s) — 



COs(E (u— Uq)) 

sin(E (u—uq)) 



and (f)(u, s) 



— sin(E(u—uo)) 

COs(E (U — Uq)) 



, and the m function 



would be the constant +£. In the general case the m function is known to be an analytic function 
mapping the upper half-plane to itself. Here, as the one-dimensional space of square-integrable 



solutions has basis 



A(u,s) 
B(u,s) 



we have m{E) 



-B aa (s) 



We had already argued from the evaluator 
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formula j8j that this meromorphic function has positive imaginary part on the half-plane $s(E) > 0. 
As our m-function is meromorphic and real-valued for E real, the spectral measure which is given 
via the formula v{a,b) = lim e _>o+ — f $s(m(E + ie))dE (under the condition v{a} + v{b} = 0) is 
a sum of Dirac measures and the spectrum is thus purely discrete, having its support on the zeros 
E n of Aa (h + iE). We obtainQ 



v{dE) 



E Bgg ( 2 
i A' 



iE„ 



iA' ao {\ + iE n ) 
From the general theory the Hilbert space L 2 ((uq, +oo) 



5(E - E n ) dE 



(22) 



C ; du) admits an orthogonal basis (T„ 



given by the two-component vectors T n (u) = <p(u 
measure formula (|22[) . norms \T n (u)\ 2 du = 
J2 n ^iKr^"" 1 dE and the Plancherel formula is 



iE n ) which have, according to the spectral 
The spectral measure is v(dE) = 



B a „(±+iE n ) 



/ \T(u)\ 2 du = Y / 



\T{E n t 



r\T n (u)\ 2 du 



T{E n ) 



T(u) ■ 4>(u, - + iE n ) du 



Let us consider the vectors 
Z n (u) = 



2A(u,±+iE n ) 
2B(u,±+iE n ) 



l«>n (u) = 2£(lt , - + iE n ) T n (u) 



(23) 



(24) 



We deduce from what has just been stated regarding the T„'s that the Z n 's are an orthogonal basis 
of L 2 ((u a , +oo) 



I 2 ; ^du) with Hilbert norms in this space: 



(Z n \Z n ) = 2iA' ao (± + iE n )B ao {\ + iE n ) 



(25) 



From Theorem [5] we know that Z n in the "JF(s)-picture" corresponds precisely to the evaluator 
Zl° +iE . Equations |2"5|) and (|10[) do indeed match. Furthermore the two-component vector Z n (u) 
solves the differential system on (mq,+oo) for the eigenvalue E n . In conclusion the differential 



operator 





-1 



d 

du 







(j,(u) 




q(u) 
f}{u) 



on (uq, oo) with the boundary condition 



acting on pairs 

a(u(j) — is self-adjoint with a purely discrete spectrum given by the imaginary parts E n of the 
zeros of A ao (ao = exp(-u ))- We notice, as ZT +lE (| + iE) = 2B aa {\ + iE n ) ^e-eT^ > tnat tnc 
eigenvectors T n {u) which are such that T n (uo) = [?] correspond exactly in the ^-"(s)-picture to the 
functions s 



A ao {±+iE) 
E—E n 



Let us denote by Hq° the self-adjoint operator given by the differential system (fl"5j) on (u , +oo) 
with boundary condition [°] at uo = log(ao). With reference to the comment after Theorem [HI we 
observe that its eigenvectors (in the J-(s) picture) s i— > (A ao (p) — 0) are not located in the 

domain of Hq (which multiplies by E), and indeed A ao (s) is not square-integrable on the critical 
line, as can be seen directly or from the presence of the Dirac at a in A (x) or the behavior of 
A ao (x) as x — > +oo. To discuss this further let us mention: 



derivatives are computed with respect to s - 



+ i 



E so iA' ao {\ + iE n ) is real. 
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Proposition 9. If f is in K a and is in the domain of Hq in L 2 , then i?o(/) is in K a . The 
restriction of Hq to K a has deficiency indices (1,1). The self-adjoint operator Hq is an extension 
of this restriction. 

Let / in K a be such that i?o(/) is hi L 2 . In the M-picture the two-component vector T repre- 
senting / must have its components absolutely continuous, and the action of Hq on T is the one 
given as a differential operator. As / is in K a , the components of T vanish identically for u < log(a) 
so their derivatives also. Hence Ho(f) £ K a . The statement is also easy to establish directly in 
the ^-picture using the formula C(x-£; + \)f = ~(x-£; + \)C(f) (in the distribution sense). Let us 
note though that by necessity the function T(u) being continuous will be such that a(log(a)) = 
and /3(log(a)) = 0. The domain of the restriction of Hq to K a , in the u-picture is thus given by the 
absolutely continuous square integrable functions, with both components vanishing at log(a), and 
whose images under the differential operator are square-integrable. This shows that the operator 
Hq is indeed an extension of this restriction. From well-known theorems |18 [ I16 j . the deficiency 
indices of such a differential operator on a semi-infinite interval, in the limit-point case at infinity, 
are (1,1) with the one-parameter family of possible boundary conditions at log(a) giving the self- 
adjoint extensions. We can confirm directly in the J-"(s)-picture the computation of the deficiency 
indices. Indeed i(Ho — ~) is multiplication by s. As sF(s) vanishes at (here T is taken in a K a , 
so is also defined away from critical line) it is perpendicular (if T is in the domain of Hq) to the 
evaluator Zq . Conversely if Q in K a is perpendicular to Z^ , the quotient \G{s) is an entire function 
and it is easy to prove that it corresponds to an element of K a (see [8] for details). So the range 
of Hq\ k — | has codimension 1, and similarly for the range of Hq\ k + |. This confirms that the 
deficiency indices are (1,1). 

6 Continuous spectrum 

Definition 10. We let L a be the perpendicular complement of K a in L 2 . In the x picture one 
has L a = L 2 (0,a;dx) + C(L 2 (0, a; dx)). In the u picture L a is the space of the square integrable 
two-component vectors supported on (— oo, log(a)). 

Again, from the M-picture we see that if / is in L a and is in the domain of Hq in L 2 then -/?o(/) is 
again in L a . And we again have a one-parameter family of self-adjoint extensions of the restriction 
of Hq to L a . But all these extensions will have continuous spectra, as the potential function yu(u) 
vanishes exponentially as u — > —oo. Summing up: 

Proposition 11. If f is in L a and is in the domain of Hq in L 2 , then i?o(/) * s ^ n L a . The 
restriction of Hq to L a has deficiency indices (1,1). Its self-adjoint extensions correspond to a 
one-parameter family of boundary conditions at log(a) and each has a purely continuous spectrum. 
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To obtain a self-adjoint operator on all of L 2 we thus need to choose the boundary condition at 
uq. A specific choice will emerge from the following procedure. The explicit spectral expansion will 
be given in the next section, after having explained how the results obtained in the present section 
fit into the framework considered by Sierra in |23j . once we go back to the a;-picture. 

We consider Hq as a differential operator in the u picture. The bound states ^"i 2 ~t lE ^ are there 



E-E n 

the vectors denoted T n (u), which have support in u > uq, and verify lim u ^ Mo T n (u) = [?] 
obey for — oo < u < +00 the differential equation, say as distributions: 



f- 

du 

-4-. 



T n (u) 



fi(u) 
fi(u) 



T n (u) = E n T n (u) 



S(u - Uq) 





We will write (recalling the use of 4r in the u picture): 



and shall convene that 
(S \T n ) = and (5i\T n ) 
geneous equation: 



£ 



S (u) 



(So\T) 
(Si\T) 
_ 1 



25(u — Uq) 



= jK. 



Si(u) 





25(u — Uq) 



They 



(26) 



(27) 



■ + lim u _ +u - )T{u) when the limits exist. So here 
This transforms the inhomogeneous equation into a linear homo- 



du 







T n (u) 



n(u) 
H(u) 



T n (u) = E n T n (u) + (5i\T n )S {u) 



(28) 



In order to obtain a formally symmetric expression we add one more term, obtaining in the end the 
following perturbation of the differential operator Hq: 



H(T) = H (T)-(S 1 \T)Sq-(Sq\T)S 1 



(29) 



According to the remark made after Theorem [SJ as a corollary to the fact that the functions 
and -^Li) are square-integrable on the critical line, when T is in the domain of Hq its components 
are continuous functions and evaluation at uq is a continuous linear form (given by equation (|18p b 
As a bilinear form, the perturbation term is thus bounded relative to (T\T) + (Hq(T)\Hq(T)). 

Let us consider for a given E the equation H(T) = ET on the full line. So T must be of the 
form T(u) — L(u)l u<Uo (u) + R(u)l u>Uo (u) where L and R verify on the full line (as differential 
equations) H (L) = EL, H (R) = ER. Writing L(u ) = [\°] and R(u ) = [%] we obtain from 
((2^|) the matching conditions: 

- lido + lo^i + riS - r S! = (li + n)5 + (lo + r )Si h = r = (30) 

We thus have for each E a two dimensional space of eigensolutions, which are of the form: 

T E (u) = XL E (u)l u<Uo (u) + /j,R E (u)l 

u>u { u ) (31) 

where L E verifies the differential equation Hq(L e ) = EL E and the condition L e (uq) = [ J] and R E 
verifies the differential equation Hq(R e ) = ER E and the condition R e (uq) = ["]. 
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As we saw in the last section the vectors Re(u)1 u>Uo (u) lead to a well-defined spectral expansion 
of K a where in fact only those among the Re which are square-integrable contribute, the suitable 
25's being the (imaginary parts of the) zeros of the entire function A ao (s). On the other hand, 
on the basis of standard theorems of spectral expansions [TBI EE E] the vectors Le(u)1 u<Uo (u) 
for E real are the generalized eigenvectors of the self-adjoint extension of the restriction of Hq to 
(— oo, log(a)), given by the boundary condition [q] at uq. 

We now extend Hg to all of L 2 defining it on L a to be this self-adjoint extension of Hq\ l • In 
the end we have thus been led by this procedure to a well-defined self-adjoint operator on L 2 whose 
spectrum has on one hand a discrete part corresponding to K a and on the other hand a continuous 
part corresponding to L a . Its eigenvectors are the solutions of the equations: 

H (T)-(8 1 \T)S -(5o\T)8 1 =ET (32) 

But only those solutions that are square-integrable as u — > +00 have to be retained, those solutions 
being: on one hand the vectors previously denoted T n (u) on (uq, +00), on the other hand a contin- 
uous spectrum of generalized elements living on (— 00, uq)- The spectral density will be determined 
in the next section. Here we now explain how the x picture looks like. 

According to Theorem [5] equation (|18p. 5q(u) corresponds in the !F{s) picture to 2.4 ao (s) and 
Si(u) corresponds to 2B ao (s) (here we are strictly on the critical line, these functions are real- valued) . 
In the f(s) picture, Sq(u) corresponds thus to 2A a (s) and S±(u) to 2B a (s). In the x picture Sq and 
£1 correspond to the distribution 2A a (x) and 2B a (x). The operator H takes the form: 

H(f) = -i(x± + \)f - (2B a \f)2A a - (2A a \f)2B a , (33) 

where the "scalar products" (2A a \f) and (2B a \f) have to be defined in a suitable manner if / does 
not belong to the domain of Hq. We will switch to another linear combination allowing an easier 
discussion and having the format considered by Sierra in [23] . 

Using E a = A a — iB a and F a = C(E a ) = A a + iB a we rewrite equation ([53"|) as 

H(f) = -i(x-£- + -)/ + 2^(^17)^ - 2i(F a \f)E a (34) 
ax 2 

The (real valued) distribution E a differs from a Dirac y /r aS a (x) by a square integrable function 
supported on [a, +00), and its Cosine transform F a = C(E a ) = A a + iB a is a function supported on 
[a, +00) which differs from y / a2cos(27raa;) by O(^) at +00. When we try, directly in the x picture, 
to solve the eigen-equation H(f) = Ef we have to give a meaning to (E a \f) where / may have 
a jump at a (we may give an imaginary part to E and look for a square-integrable / to simplify 
the discussion). The correct choice is to take (S a \f) for those / to be the average between the left 
and right values at x = a. And (F a \f) needs to be defined suitably as well. The discussion in the 
u-picture was technically easier and we shall defer to another publication the more detailed direct 
treatment of (|34jl in the x-picture. 
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Let us mention that if we repeat the whole procedure but starting from the self-adjoint operator 
on K a corresponding to the zeros of the B a function we end up with the following operator 

H(f) = -i(x-f + hf + 2i{F a \f)E a - 2i(E a \f)F a (35) 
dx 2 

which differs from the previous one only by the sign of the perturbation. 
The general equation considered by Sierra in [23l IV] is 

H(f) = -i{x± + \)f + i(HfWi - (36) 

Thus (|34p and (f3"S"|) are indeed special cases. Furthermore Sierra considers the "duality" condition 
that and ip2 are a Cosine (or Sine) transform pair. Our special instances do have this property. 
Finally in [23l VI] the case ip2(x) = V2^/a S a (x) and tpi = C(ip2) is devoted special attention. 
Whether ip2(x) = V2yfa 5 a (x) leads to bound states is not clear. 



7 More on the scattering spectrum 



We want to provide more information on the "scattering" part of the spectrum. When E is not real 
the differential system (fT5j) has (up to a multiple) exactly one solution which is square-integrable 
at — oo. Fortunately, this solution is known explicitely. It was determined in [7], as well as the way 
BOu's] f° r E rea ^ ^ S ex P resse d as a linear combination of the boundary values of such solutions 
and their complex conjugates. 



a E {u) 
Pe{u) 



of the differential system which 



Theorem 12 ((?]; Theorem 16 of [HI])- The eigensolution 

is square-integrable at —oo for ^(E) > (Ui(s) < lj is provided by the following formulas (where 
a ~ exp(it) ): 



otE(u) — J(u, — + iE) — t 



-iEv 



4>a( x ) x * d x 



Pe{u) = K(u, -+iE) = ie~ ltjU + ie? u / 4>~ (x)ar s dx 
2 Jo 



The following identity holds: 



2A(u, \ + iE) 
2B(u, | + iE) 



= 7r-*r(-i 



J(u, h + iE)' 
K(u±+iE) 



V 2 ; 



J(u, \-iE) ' 
-K(u,\-iE) 



(37) 
(38) 

(39) 



We will also use the notation j a {s) = J(u, s) and k a (s) — K(u, s), for u — log(a): 



J(u,s) — j a (s) = a 2 3 — ya / (f)J(x)x s dx 

Jo 

^ i [ a 

K(u, s) — k a (s) — ia 2 ~ s + i\fa I <fi~ (x)x~ s dx 



(40a) 
(40b) 
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The functions <j)J (x) and 4>~ (x) being entire and even, we see that j a and k a are meromorphic in 
the entire complex plane with possible poles at s = 1, s = 3, s = 5, ... Of course they still provide 
after analytic continuation a solution of the differential system with respect to log(a). We have 

2A a (s) = 7T-*r(!)£( fl ) + 7^1X^)^(1 - S ) (41) 

As ^la(s) is entire, if j a does not have a pole at 1 + 2m (m G N), necessarily the second term on 
the right-hand side is forbidden to have a pole there, which means that necessarily one must have 
a zero of j a at —2m. We will now argue that neither j a nor k a can have zeros in the half-plane 
5ft(s) < | (which corresponds to 3(i?) > 0). This will prove that j a (and similarly k a ) does have 
poles at s = 1, s = 3, s = 5, . . . 

Indeed as in the previous section we apply Hermann Weyl's theory ([UJ §9], [16j §3]), but on the 
interval (™oo, uq). We maintain the same notations ip(u, s) with boundary condition [J] and (/'(u, s) 
with boundary condition [J] but convene that the m-function for ^s(E) > is determined by the 
condition that —m(E) , ip(u, s) + (/>(it, s) is square integrable at — oo. Thus, we have the identity: 

m > =* m(E) = =g^f = ^ (42) 

It is a general theorem that the m function is an analytic function mapping the upper half-plane 

to itself. So if k ao (s) has a zero for 5R(s) < i, then j ao (s) has to vanish also for the same sq. But 

KOa^o) * s s °l u ti° n °f a first order differential equation with respect to u. If it vanishes for some 

u it has to vanish for all it's. From (|40ap . (|40b[) however it is easily argued that j a (s) ~ a^~ s as 

a -> 0+, and j£(s) ~ Q _ «a^ s (3?(s) < 1)@ 

.4(u, s) J(u, s) 
B(u,s) K(u,s) 

perhaps first for < 3?(s) < |, using the asymptotic of J and K as a — > + and the identity 
This gives: 

W {s)=i^T{^-) (43) 

Analytic continuation removes the restriction on s. This confirms that J or K (in fact both as we 
saw) must have poles at 1, 3, . . . 



The Wronskian W(s) 



is independent of u. We can evaluate it at u 



The Wronskian W\(u,s) 
equation 



J(u, s) J(u, s) 



K(u, s) K(u, s) 



Q(— J(w, s)K(u, s)) verifies the differential 



4-Wi{u,s) = S(E) (|J( U ,s)| 2 + \K(u,s)\ 2 ) (44) 
au 

For Sft(s) < 1, we have Wi(u,s) ~ o _v0 a 1_2!R ( s ). For 3?(s) < |, we have 3(E) > and deduce from 
((SI) that S$(—J(u,s)K(u,s)) > which confirms the m-property of the ratio ~k(™'s) ( an d the fact 



8 the equations ]40al l Il40bl l also give easily for a given a asymptotic expansions of j a (s), k a (s) in powers of — in 
the critical strip, and this gives thus corresponding information on A a , B a and E a . 
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that none of J and K can vanish in this half-plane, whatever the value of u). Additionally, for 
3?(s) = \ we obtain: 



Proposition 13. For s = | + iE, figl and — oo < it < +00 i/iere holds 

<s(-J{u,s)K(u, s)) = 1 
Hence the functions J(u,s) and K(u,s) do not vanish on the line 5R(s) = i. 



(45) 



value — 2i of the Wronskian 



from the behavior at —00. The spec- 

_ -./(mq.s) 



Alternatively, we could also have deduced (not only for E real but for all E £ C) the constant 

J(u±+iE) J{u,\-iE) 
K(u,\+iE) -K(u,\-iE) 

tral measure v{a,b) = lim e ^ + ^ f a $l(m(E + ie))dE, m(E) = -g^^y is thus purely absolutely 

continuous: 

<dE) = ^ (46) 

The spectral expansion takes the following form: we send T{u) 



f3(u) 



to 



r«o 2. 1 

T{E)= I a(u)i/;o(u,-+iE) + p(u)ip 1 (u,- + iE)du 
-00 ^ ^ 



(47) 



where ip(u,s) — is the unique solution of the differential system verifying the boundary 

condition [J] at uo,lj and this transform is isometric with the following Plancherel identity: 



U0 |a( M )| 2 + |/?( M )| 2 $ = i- / |T(£)| 2 



2tt 



dE 



\K(u ±+iEW 



(48) 



The two linearly independent eigensolutions of the differential system are 
with Wronskian — 2i. Hence: 



J(u, \ - iE) 



J{u±+iE) 
K(u 7 \ + iE) 



and 



-K(u, ± - iE) 



-i K(uq, \ — iE) 



J(u, h+iE)' 
K{u,\+iE) 



-iK{uo,\ + iE) 



J{u, \ - iE) 
-K{u,\-iE) 



(49) 



1 



1 



Mu, r + iE) = 9f K> , = - i£)J(u, - + iE) 



^{u, - + = 3 if(u , - - i£)A(n, - + iE) 



(50) 
(51) 



9 as s = = + iE, E £ R, this solution is real valued. 
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